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Chapter  1 
INTRODUCTION 

The  release  of  volatile  organic  compounds  (VOC's)into 
the  atmosphere  has  come  under  fire  by  the  Environmental 
Protection  Agency  (EPA)  and  state  regulatory  agencies  in  the 
past  several  years.   These  releases  come  from  a  wide  variety 
of  sources  from  gasoline  filling  stations  to  industrial  and 
military  operations.   There  are  also  countless  natural 
sources  of  VOC's  including  volcanos  and  forest  fires. 
Environmental  groups  have  had  a  great  influence  on 
government  agencies  resulting  in  the  enactment  of  numerous 
statutes  and  regulations  to  reduce  these  releases  over  an 
approximate  five  (5)  year  period. 

Enforcement  of  environmental  regulations  is  carried  out 
by  the  individual  states  which  impose  their  own  specific 
regulations  which  in  some  cases  are  even  more  stringent  than 
those  of  the  federal  government.   A  case  in  point  is  the 
state  of  California  where  regulations  are  resulting  in  rapid 
actions  on  the  part  of  corporations  to  comply  or  face  harsh 
penalties.   These  actions  include  both  a  switch  to 
alternative  manufacturing  materials  and  employment  of  new 
pollution  control  devices. 

A  new  pollution  control  technology  currently  under 
development  utilizes  ultraviolet  (UV)  light  to  break  down 
the  volatile  organic  compounds  (VOC's)  contained  in  an 


2 

airstream.   One  pollution  control  equipment  manufacturer  has 
already  employed  this  method  in  an  installation  in 
California  (see  Figure  (1) ) .   The  process  involves  injecting 
ozone  into  a  mixture  of  gaseous  VOC  s  which  have  been 
radiated  by  UV  light.   In  the  presence  of  oxidants  and  UV 
light,  the  air/VOC  mixture  undergoes  a  series  of  photo- 
dissociation  and  radical  oxidation  reactions  which  destroy 
the  VOC's.   Despite  this  pollution  control  equipment's 
success  in  controlling  a  wide  range  of  airborne  VOC's,  the 
extent  of  VOC  reduction  due  to  photolysis  in  the  photolytic 
reactor  section  of  the  system  is  unknown.   Our  lack  of 
information  on  the  effectiveness  of  this  photolytic  reactor 
as  shown  in  Figure  (2)  is  due  to  proprietary  information 
retained  by  the  manufacturer,  insufficient  test  data,  and 
lack  of  an  adequate  model  to  predict  its  performance. 
There  have  been  some  analytical  and  experimental 
studies  reported  in  the  literature  on  photochemical  reactors 
(ie.,  Y.  Harano  and  T.  Matsura  (1972)).   However,  these 
studies  are  generally  directed  towards  simplified  geometries 
which  can  be  modeled  analytically  (ie.,  axial,  elliptical, 
and  columnar  reactors)  and  which  assume  well-mixed 
conditions,  low  light  absorption,  uniform  flow,  and 
negligible  diffusion.   Experimentation  has  been  directed 
towards  determining  the  deviation  between  these  highly 
simplified  models  and  actual  reactor  performance  in  order  to 
determine  the  effects  of  factors  omitted.   Few  if  any 
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Figure  (2)  -  VOC  Pollution  Control  Device  Photolytic  Reactor 


5 
studies  can  be  found  in  the  literature  which  attempt  to 
quantitatively  model  photochemical  reactors  of  complex 
geometries  and  flow  fields  using  numerical  methods.   This 
report  provides  a  numerical  model  which  can  be  used  to 
quantitatively  describe  the  photo-dissociation/radical 
oxidation  of  VOC's  in  the  photolytic  reactor  of  Figure  (2). 

A  secondary  but  equally  important  objective  for  this 
analysis  is  to  develop  an  analytical  tool  which  can  be  used 
to  optimize  photo-reactor  geometry.   This  numerical  model, 
once  verified  with  experimental  data  (in  progress  at  Penn 
State)  and  expanded  to  incorporate  additional  reactants,  can 
be  utilized  to  predict  optimum  geometry  for  the  photolytic 
reactor  section  of  a  pollution  control  device. 

Finally,  this  model  will  assist  the  designer  in 
determining  when  the  use  of  a  photolytic  reactor  in  a 
pollution  control  device  is  appropriate.   Since  there  are  a 
multitude  of  VOC's  which  vary  widely  in  their  ability  to 
react  when  exposed  to  UV  light,  cost  effective 
implementation  of  a  photolytic  reactor  may  be  highly 
dependant  on  the  specific  VOC's  expected  in  an  airstream. 
For  example,  the  photolytic  reactor  may  not  provide  any 
additional  benefit  in  a  pollution  control  device  which 
treats  specific  emissions  from  a  military  or  industrial 
painting  operation.   This  numerical  model  may  indicate  when 
a  significant  cost  savings  is  possible  by  omitting  the 
photolytic  reactor  section  and  could  also  predict  when  the 


6 
photolytic  reactor  can  function  as  the  sole  VOC  reduction 
device  in  a  pollution  control  system. 

1 . 1   Literature  Review 

A  review  of  literature  on  the  subject  of  design  and 
analysis  of  photochemical  reactors  shows  that  very  little 
work  has  been  accomplished  in  the  area  of  numerical 
modeling. 

Jacob,  S.  M.  and  Dranoff,  J.  S.  (1969)  conducted 
experiments  on  light  distribution  within  an  elliptical 
photoreactor  and  compared  these  results  with  some  closed 
form  solutions.   They  found  that  the  closed  form  solutions 
are  unable  to  accurately  predict  light  intensity  within  an 
elliptical  reactor.   Additional  studies  by  Jacob,  S.  M.  and 
Dranoff,  J.  S.  (1970)  concluded  that  numerical  means  are 
required  to  accurately  determine  the  light  intensity 
distribution  within  a  well-mixed  elliptical  photoreactor. 
While  pointing  out  the  necessity  to  use  numerical  methods  in 
the  analysis  of  photolytic  reactors,  neither  of  the  above 
studies  attempted  to  deal  with  reactors  containing  non- 
uniform flow  fields  and  spatially  varying  reactant 
concentrations.   In  both  studies  only  the  case  of  simple 
elliptical  geometry  with  laminar  flow  and  well-mixed 
conditions  was  considered. 

Work  by  Ragonese,  F.  P.  and  Williams,  J.  A.  (1971)  on 
the  scale-up  of  a  laboratory  photoreactor  model  to  meet 
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actual  chemical  production  requirements  outlines  the 
difficulties  in  utilizing  both  dimensional  analysis  and 
empirical  data.   Although  success  has  been  achieved  in  this 
scale-up  process,  the  methods  used  apply  to  simplified 
geometries  with  laminar  flow  and  well-mixed  conditions. 
This  study  points  out  the  need  for  numerical  modeling  to 
enable  a  designer  to  scale  photoreactors  to  meet  production 
needs  or  to  achieve  sufficient  pollutant  reductions  in  a 
pollution  control  device. 

The  literature  review  provided  by  Harano,  Y.,  and 
Matsura,  T.  (1972)  describes  the  many  difficulties  involved 
in  the  design  of  photochemical  reactors.   They  summarized 
that  even  for  the  simplified  geometries  which  have  been 
studied  extensively  to  date  including  cylindrical,  parallel 
plate,  annular  and  elliptical  reactors,  the  actual 
concentration  gradients,  light  intensity  distributions  and 
flow  conditions  greatly  influence  overall  reaction  rates 
within  the  photoreactor . 

All  the  literature  reviewed  on  the  subject  of 
photolytic  reactor  design  points  to  the  need  for  numerical 
reactor  modeling  which  incorporates  both  photochemical  and 
fluid  mechanics  principles  and  the  latest  numerical  solution 
techniques . 


Chapter  2 
PRELIMINARY  MODELING  ASSUMPTIONS 

Due  to  the  great  difficulty  in  modeling  the  entire 
reactor  under  investigation  including  all  possible  chemical 
reactions,  this  model  describes  only  the  photo- 
dissociation/radical  oxidation  (photolysis)  of  formaldehyde 
in  a  rectangular  flow  field  perpendicular  to  the  axis  of  a 
single  UV  lamp  as  shown  in  Figure  (2) .   Formaldehyde  has 
been  chosen  for  this  model  since  its  photochemistry  has  been 
studied,  documented  and  verified.   Although  this  model 
applies  only  to  the  photolysis  of  formaldehyde,  it  is 
general  enough  to  incorporate  the  chemistry  of  other 
airborne  components  as  well  as  combinations  of  other 
photochemically  reactive  solvents  in  the  future.   In 
addition,  the  principles  required  to  expand  this  numerical 
analysis  to  additional  materials  and/or  chemical  reactions 
have  been  included  in  this  report. 

The  following  assumptions  have  been  made  to  simplify 
the  analysis: 

A.  Constant  fluid  properties  (ie. , density, 
temperature,  diffusion  coefficient) . 

B.  Negligible  thermal  effects  (ie. , solution  of  energy 
conservation  equations  not  required) . 

C.  Two-dimensional  flow  field. 
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D.  Weak  absorbance  of  UV  light  (ie.,  light  intensity 
is  independent  of  formaldehyde  concentration) . 

E.  Photo-dissociation/radical  oxidation  of 
formaldehyde  only.   No  other  reactions  between  formaldehyde 
and  other  compounds  are  included  except  those  required  to 
complete  the  photo-dissociation/radical  oxidation  of 
formaldehyde. 

F.  Steady,  time-averaged  flow  field  (ie,  no 
oscillation  of  vortex  shedding  in  the  wake  region) . 

G.  No  wall  losses  (ie.,  chemical  reactions  which 
sometimes  occur  due  to  collisions  of  molecules  into  the 
reactor  walls  or  the  UV  bulb  itself) . 

H.   No  reflection  of  UV  light  from  the  reactor  walls. 

The  above  assumptions  are  considered  to  be  reasonable 
ones  for  use  in  this  initial  numerical  model.   In  the 
future,  factors  such  as  wall  losses  and  reflection  can  be 
introduced  into  the  numerical  simulation  to  increase  its 
accuracy. 
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Chapter  3 
MASS  CONSERVATION  WITHIN  THE  PHOTOLYTIC  REACTOR 

Given  the  assumptions  provided  in  Chapter  2,  the 
differential  equation  which  describes  material  concentration 
in  the  photolytic  reactor  can  be  written: 

-v°Vc  +  DV2C  +  Q.     =   dC/dt  [1] 

where,  C  =  concentration  of  formaldehyde  (PPM) . 

v  =  flow  velocity  (with  components  vx,v  ) (m/s) . 

D  =  formaldehyde  diffusion  coefficient  (m/s) . 
£1   =   generation  of  formaldehyde  (PPM/s) . 

The  first  element  of  Equation  [1],  the  convection  term, 
involves  the  fluid  mechanics  of  the  reactor.   It  combines 
concentration  changes  along  a  given  direction  with  the  fluid 
velocity  to  result  in  a  time  rate  of  change  in 
concentration . 

The  second  term  in  the  above  equation  is  the  diffusion 
term.   It  describes  the  time  rate  of  change  of  concentration 
entering/leaving  the  control  volume  due  to  concentration 
gradients  at  the  element  surface.   Note  that  this  term  is 
dependant  on  the  second  spatial  derivative  of  C  multiplied 
by  a  diffusion  constant  and  is  relatively  small  with  respect 
to  the  velocities  in  the  flow.   To  achieve  the  highest 
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accuracy  possible  for  all  possible  concentration  gradients 
and  flow  velocities,  this  term  is  included  in  the  model. 

The  third  term  in  this  equation  describes  concentration 
changes  due  to  the  photolytic  action  of  the  UV  bulb  and 
other  chemical  reactions.   For  this  model  the  "generation 
term"  simplifies  to  the  following: 

CI  =   -KC  [2] 

where,  K  is  the  overall  specific  absorption  rate  constant 
(overall  rate  constant  for  the  breakdown  of  formaldehyde  by 
UV  light) . 

The  final  term  is  the  resultant  time  rate  of  change  for 
the  element  which,  when  included  yields  a  transient 
solution.   Since  the  transient  solution  has  practical  value 
both  in  predicting  the  start-up  concentrations  within  the 
reactor  as  well  as  in  examining  time  varying  inputs,  the 
time  dependant  term  has  been  included.   In  addition  to  the 
above  benefits,  a  transient  solution  must  be  found  in  order 
to  obtain  convergence  for  the  numerical  solution  technique 
used.   This  is  required  since  it  is  difficult  for  the 
iterative  numerical  technique  (Gauss-Seidel)  to  jump 
directly  from  zero  time  to  an  infinite  solution  (steady- 
state  solution)  without  some  intermediate  steps  in  time. 
For  the  case  of  a  multi-component  flow  field  which 
includes  chemical  reactions  as  well  as  photolytic  reactions, 
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Equations  [1]  and  [2]  can  be  written: 

-v°VCi  +  DiV'Ci  +  Qi  =  dCj/dt  [3] 

Q^  =  -K.j_C.j_  +  other  reactions  (see  Appendix  B)      [4] 

Analysis  of  the  multi-component  case  requires  the  solution 
of  n  simultaneous  Equations  ([3]  and  [4])  corresponding  to 
the  total  number  of  reactants.   Note  that  Appendix  B 
contains  numerous  additional  equations  involving  chemical 
and  photochemical  reactions  between  many  possible  airborne 
constituents.   In  order  to  simplify  this  analysis  only  the 
photolysis  of  formaldehyde  will  be  considered  (ie.,  no 
reactions  between  formaldehyde  and  other  airborne 
constituents  will  be  considered) . 

In  the  next  chapter  the  determination  of  velocities 
required  in  the  convection  term  of  Equation  [1]  is 
discussed. 
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Chapter  4 
FLUID  MECHANICS  PRINCIPLES 

For  the  case  of  turbulent,  viscous  flow  around  the 
spherical  bulbs  shown  in  Figure  (2) ,  a  complete  solution  of 
the  Navier-Stokes  equations  is  required  to  obtain  accurate 
fluid  velocities.   Due  to  the  difficulty  in  solving  the 
complete  Navier-Stokes  equations,  a  commercially  available 
computational  fluid  dynamics  program  has  been  employed  to 
determine  the  velocities  in  a  rectangular  region  around  a 
single  bulb.   Figure (3)  shows  a  standard  K-E  method  solution 
of  the  Navier-Stokes  equations  for  flow  of  air  around  a 
circular  cylinder  generated  by  the  FLUENT  program.   The 
shape  of  a  0.03m  diameter  circular  cylinder  in  a  0.05m  x 
0.20m  flow  field  has  been  approximated  by  a  series  of  grid 
points  (100  wide  by  25  high) .   The  program  calculates  x  and 
y  velocities  at  every  grid  point  in  the  flow  for  input  into 
the  convection  terms  of  Equation  [1] . 

The  following  parameters  were  input  to  the  FLUENT 
program  (properties  of  air/solution  parameters) : 

Properties  of  air  only  assumed  (no  air/form,  mixture) 

Length  of  domain  (x  direction,  m) -  0.2 

Height  of  domain  (y  direction,  m) -  0.05 

Diameter  of  Bulb  (m) -  0.03 

Grid  size  (x  dir,  y  dir)-  100  x  25 
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Reynolds  Number,  ReD  -  4,500 
Inlet  Temperature  (°C)-  35 
Inlet  Pressure  (Kpa)-  101.33 
Inlet  Density  (kg/m3)-  1.293 
Inlet  Viscosity  (kg/m-s)-  1.72  x  10~5 
Inlet  Turbulence  Intensity  (%)-  10.0 
Inlet  Velocity  (m/s)-  2.0 

Since  it  is  beyond  the  scope  of  this  report  to  examine 
the  detailed  accuracy  of  the  FLUENT  program  output,  these 
results  will  be  assumed  to  be  sufficiently  accurate  for  the 
purpose  of  making  design  decisions.   There  are  a  variety  of 
models  and  techniques  which  can  be  employed  to  increase  and 
verify  the  accuracy  of  the  flow  analysis.   However,  since  it 
is  the  purpose  of  this  report  to  model  the  photolytic 
reactor  in  a  fashion  which  well-exceeds  the  accuracy  of 
simple  plug  flow  models  which  can  be  found  in  the 
literature,  the  unmodified  FLUENT  output  is  used. 
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Chapter  5 
CHEMICAL  PRINCIPLES 

The  chemistry  associated  with  the  photolytic  reactor 
under  consideration  here  is  much  the  same  as  that  of 
atmospheric  chemistry  with  the  exception  that  energy  is 
added  with  a  mercury  lamp  which  emits  several  discrete 
wavelengths  (predominantly  253nm)  versus  sunlight  which 
contains  these  and  other  wavelengths.   The  chemistry  of  the 
photolytic  reactor  can  be  broken  down  into  two  categories 
including  photochemistry  and  free  radical  chemistry. 

5 . 1   PHOTOCHEMISTRY 

Photochemical  processes  initiated  by  UV  radiation  begin 

with  the  absorption  of  a  photon  by  a  molecule  as  follows: 

K 
A  +  hi)  ->  A*  [  5  ] 

where,  A*  is  an  exited  A  molecule  and  hi)  is  a  photon  of  UV 
energy.   Upon  excitation  the  A*  molecule  can  undergo  the 
following  reactions: 


Kl 
Dissociation:      A*  — »  B  +  C  [6] 

(A*  broken  down  into  stable  molecules  by  itself) 

K2 
Direct  Reaction:   A*+B  — »D+E  [7] 

(A*  combines  with  B  to  form  two  stable  molecules) 
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K3 
Fluorescence:      A*  — >A  +  hv  [8] 

(A*  returns  to  orig.  state  after  release  of  energy) 

Collisional  K4 

Deactivation:       A*  +  B  — >A  +  B  +  hD  [9] 

(A*  returns  to  orig.  state  after  collision  with  B) 


Each  one  of  the  above  reactions  has  a  probability  of 
occurrence  called  a  quantum  yield,  <J>j_ .   For  example,  a 
quantum  yield  of  4^  =  0.5  implies  that  there  is  50%  chance 
that  A*  molecules  will  undergo  dissociation  (Equation  [6]) 
while  there  is  a  50%  chance  A*  molecules  will  directly 
react,  undergo  fluorescence  or  deactivate  from  a  collision 
Since  all  the  A*  molecules  must  undergo  one  of  the  above 
reactions : 


l 

hi 


E  <&±  =  i  [io 


The  production  of  A*  molecules  is  given  by  the 
following  rate  equation: 


dTA*1  =  K[A]  [11] 

dt 


In  Equation  [5],  since  one  A*  molecule  consumes  one 
molecule  of  A: 


dTAI  =  -K[A]  [12] 

dt 
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The  symbol  K  is  the  overall  "specific  absorption  rate" 
constant  (also  given  in  Equation  [2])  and  [A*],  [A]  are  the 
concentrations  of  A*  and  A,  respectively.   Note  that  the 
overall  rate  constant  corresponds  to  all  possible  quantum 
yields  for  A*.   By  combining  equations  [5]  through  [9],  the 
rate  equations  corresponding  to  the  dissociation,  direct 
reaction,  fluorescence,  and  collisional  deactivation 
reactions  (termed  the  first  order  rate  equations)  can  be 
written  as  follows: 


dfAl  =  -<E>1K[A]  =  -KX[A]  (reduction  in  A  due  to     [13] 
dt  dissociation) 

d[Al  =  -^2KtA^  =  ~K2 t AJ  (reduction  in  A  due  to     [14] 
dt  direct  reaction) 

dTA]  =  -<X>3K[A]  =  -K3[A]  (reduction  in  A  due  to     [15] 
dt  fluorescence) 

dTA]  =  -d>4K[A]  =  -K4[A]  (reduction  in  A  due  to     [16] 
dt  collisional  deactivation) 


The  quantities  K-,_  (^-lK)  ,  K2  K3  and  K4  in  Equations  [13] 
through  [16]  are  termed  the  first  order  rate  constants. 
Note  that  multiplying  both  sides  of  Equation  [10]  by  K 
results  in  Kx  +  K2  +  K3  +  K4  =  K.   Therefore,  the  overall 
reaction  rate  for  the  reduction  of  A  is  equal  to  the  sum  of 
the  individual  rates  for  reactions  which  can  occur  upon 
excitation  of  A  molecules. 

For  the  general  case  of  atmospheric  chemistry,  the 
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first  order  rate  constants  can  be  determined  through  the 
following  equation: 

K-l  =  \  C(A.,T)<frL(A.,T)  I(X)dX  [17] 

G(A,,T)  =  absorption  cross  section  of  the  A  molecule  at 
operating  temp.,  pressure  and  UV  wavelength 
(cm2) . 

<!>-]_  (A.,  T)  =  quantum  yield  or  probability  that  molecule  A  will 
dissociate  upon  excitation  (recalling  that  one 
molecule  of  A  produces  one  molecule  of  A*)  (K2 
would  correspond  to  direct  reaction  rate) . 

I(A,)dA.  =  irradiance  of  incident  UV  light  (photons/cm2-sec)  . 

The  above  first  order  rate  constant,  K-j_  corresponds  to 
the  rate  at  which  molecule  A  will  undergo  dissociation  (one 
of  the  four  reactions  mentioned  above)  when  exposed  to  a 
broad  spectrum  of  light  energy  in  the  atmosphere. 

In  order  to  calculate  the  first  order  rate  constants 
for  combinations  of  specific,  discrete  wavelengths  the 
following  approximation  is  utilized: 


%  =  E  a(\±fr£)&1(\±,T)H\i)&k  [18] 

is/ 
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The  quantities  G  and  <£-]_  above  are  average  values  at  X^.      The 
quantity  I(X.i)AX,  is  sometimes  expressed  as  J(A,i)  (actinic 
flux)  which  results  in  the  following: 


n 

K 


1  =   E  a(X1,T)*L(X1/T)J(Xi)  [19] 


<=< 


5.1.1   DETERMINATION  OF  FORMALDEHYDE  PHOTOLYTIC  RATE 
CONSTANTS 

Formaldehyde  participates  in  the  following  photolytic 
reactions  (A.<300nm)  : 


Kl 
HCHO  +  hv  — >H  +  HCO  (HCHO  is  formaldehyde)         [20] 


K2 


H2  +  CO  [21] 


The  wavelength  of  UV  light  utilized  in  this  device  is 
predominately  X  =  253  nm.  The  bulb  dimensions  are  as 
follows : 

D  (diameter)  =  0.033  m 
L  (length)  =   1.5536  m 

ptotal  (total  bulb  power  requirement)  =  40  W  (J/s) 
^visible  (bulb  output  in  visible  spectrum  per 
manufacturer)  =  40%  Ptotai 

P253  (bulb  output  at  253  nm)  =  95%  Pvisible 
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As  an  example,  the  first  order  rate  constant  at  A.=253nm 
(corresponding  to  the  reaction  of  equation  [20]  above)  is 
determined  as  follows: 

E  =  he  A  =  6.625  6xl0"34J-s/photon(2.9  97  9xl08m/s) 

253x10  *m 

=  7.8509  x  10"19J/photon 

ABULB  =  7cDL  =  ft (33mm)  (1553.6mm)  =  1611  cm2 
p253  =  40W(40%)  (95%)=  14.44  J/s 

J(X=253nm)  =  _J?253 =   14.44  J/s . 

ABULB(E)   1611  cm2 (7.8509  x  10~19J/photon) 
=  1.1419  x  1016  photons/cm2s 

<T(253nm,25C)  =  3.79  x  10~21cm2/molecule  (R.  Atkinson, 

et  al.  (1992) ) 

O-l  (253nm,25C)  =0.3  (R.  Atkinson,  et  al .  (1992)) 

K±   =  a(X,T)<b1(\tT)  J(X) 

=  3.79xl0"21cm2/molecule(0.3) (1 . 1419xl016phot . /cm2s) 
=  1.30  x  10~5  (1/s) 

K  =  overall  specific  absorption  rate  (yields  both 
reactions  1,2) 
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K  =  %  +  K2  .  o(X,T)®1(X,T)  J  (A.)  +  a(A,,  T)0>2(A,,T)  j(A.) 

=  a(A.,T)  {Ql(X,?)-*&2(\rT))J(X) 
02(253nm/25C)  =0.5  (R.  Atkinson,  et  al .  (1992)) 

K  =  (3.79xl0~21cm2/molecule) (0.3+0.5)  x 

(1.1419xl016photons/cm2s) 

=  3.47  x  10~5  (1/s)  (yields  both  reactions  1,2 

(Eqns  [20],  [21])) 

5.1.2   LIGHT  INTENSITY  DISTRIBUTION  THROUGHOUT  THE  REACTOR 

The  rate  calculations  above  are  based  on  the  intensity 
at  the  bulb  surface.   This  formulation  is  based  on  the  Beer- 
Lambert  law  applied  to  weak  absorbance  conditions  such  as 
the  atmosphere  (Finlayson-Pitts,  B.  J.,  and  Pitts,  J.  N., 
Jr.  (198  6) )  .   Since  the  energy  added  via  the  bulb  per  unit 
area  and  corresponding  rate  constants  vary  indirectly  with 
the  radial  distance  from  the  bulb  surface,  Equation  [2]  must 
be  modified  as  follows: 

Q  =  -K(R/r)C  [22] 

where,       r  =  radial  distance  from  bulb  centerline. 
R  =  radius  of  bulb. 

In  addition,  the  light  intensity  in  an  actual 
photochemical  reactor  is  also  attenuated  according  to  the 
well  known  Beer-Lambert  law  (Finlayson-Pitts,  B.  J.,  and 
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Pitts,  J.  N.,  Jr.  (1986))which  is  written  as  follows: 

I  =  Io  exp{-(JNP}  [23] 

where,  Io  =  monochromatic  source  intensity (photons/cm2-sec) . 
o  =  absorption  cross  section  (cm2/molecule) . 
N  =  molecules/cm3  (molar  concentration) . 

P  =  path  length  from  the  source  (cm) . 

For  low  concentrations  and/or  low  absorbance,  as  is  the 
case  with  this  model,  the  effects  of  concentration  on  light 
intensity  can  be  shown  to  be  negligible.   For  this  reactor: 

a   =  3.79  x  10~21cm2/molecule 

N  =  1.229  x  IO15  molecules/cm3  (equivalent  to  50  PPM) 

P  =  approx.  20  cm 


Substituting  into  Equation  [23]  yields: 

I  =  Io  expi-(3.79  x  10-21cm2/molecule)  x 

(1.229  x  IO15  molecules/cm3) (20  cm)) 
=  0.99991  Io 


In  view  of  the  above  result  it  is  clear  that  the  weak 
absorbance  assumption  is  a  reasonable  one. 
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In  order  to  accurately  model  a  multi-component  air/VOC 
mixture  where  light  absorbance  must  be  considered,  the 
following  multi-reactant  Beer-Lambert  law  is  utilized: 

I  =  io  expf-ta-jN-L  +  a2N2  +  ]L}  [24] 

Using  this  relationship  an  intensity  function,  I  =  f (x,y) 
can  be  determined  using  the  actual  concentration  of 
materials  in  the  reactor  which  are  capable  of  absorbing 
light. 

The  specific  photochemistry  of  reactants  which  could  be 
considered  in  this  model  are  contained  in  Appendix  B. 
However,  for  the  purpose  of  this  initial  model  only  the 
photolysis  of  formaldehyde  is  under  consideration. 
Therefore,  the  intensity  will  be  considered  to  be  a  function 
of  the  radial  distance  from  the  bulb  surface  only  due  to  the 
weak  absorbance  of  formaldehyde  which  was  demonstrated 
above. 

5.2   FREE  RADICAL  CHEMISTRY 

As  can  be  seen  in  Appendix  B  there  are  a  multitude  of 
reactions  which  can  be  included  in  this  model  as  there  are 
numerous  materials  normally  contained  in  an  airstream  which 
are  capable  of  undergoing  photolysis  and  reacting  with 
formaldehyde.   However,  due  to  the  great  complexity  of 
including  all  these  reactions,  this  initial  model  will  be 
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only  concerned  with  the  removal  of  formaldehyde  through 
photolysis  alone.   Therefore  the  "generation"  term  of 
Equation  [1]  will  be  that  given  by  Equation  [20]  with  no 
other  reactions  included. 

Note  that  in  making  the  above  simplification,  reactions 
with  free  radicals  have  been  ignored.   Therefore,  the 
results  of  this  investigation  apply  only  to  photolytic 
reactions.   This  analysis  is  a  logical  first  step  toward  a 
more  comprehensive  analysis  that  will  include  free  radical 
reactions.   This  analysis  also  provides  valuable  insight 
into  the  photolysis  of  formaldehyde  in  a  non-uniform  flow 
field. 
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Chapter  6 
CLOSED  FORM  SOLUTION  FOR  SIMPLIFIED  FLOW  FIELD 

In  order  to  get  a  general  idea  of  the  size  of  the 
actual  pollution  control  device  required  to  reduce  a  given 
formaldehyde  input  concentration  by  a  factor  of  two  through 
photolysis,  a  simplified  flow  field  model  can  be  employed. 
Restating  the  mass  conservation  equation: 

-v°VC  +  DV2C  +  Q  =  dC/dt  [25] 

The  following  simplifying  assumptions  will  be  made: 

A.  Two-dimensional,  inviscid,  steady  flow. 

(v  =  velocity  in  x  direction  =  constant, 
velocity  in  y  direction  =  0) 

B.  Constant  radiance  intensity  throughout  the  reactor. 

C.  Formaldehyde  photolysis  only. 

D.  Diffusion  negligible. 

With  these  assumptions,  the  above  equation  simplifies  to: 

-v°VC  +   Q  =    0  [26] 

Since  the  only  reaction  involves  the  photolysis  of 
formaldehyde,  substitution  of  Equation  [2]  into  Equation 
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[26],    performing  the  dot  product  and  simplifying  yields: 


v  dC  =  -KC  [27] 

dx 


Rearranging  terms  and  solving  for  C(x)  yields 


ln(C)-ln(CQ)  =  -Kx/v  [28 

ln(C/CQ)  =  -Kx/v  [29 

C/C0  =  exp{-Kx/v}  [30; 


The  values  assumed  are 


C/CQ  =0.5 

K  =  3.47  x  10"5  (1/s)  (as  calculated) 


v  =  5  ft/s 

Substituting  in  values  and  solving  for  x  yields: 

X50  =  99,811    ft  =  x  corresponding  to  a  50%  C  reduction 
=  18. 9  miles 

This  result  is  a  very  pessimistic  one  which  predicts 
poor  performance  for  a  reactor  which  relies  on  the 
photolysis  of  formaldehyde  alone.   However,  in  examining  the 
general  solution  of  the  analytical  model  for  different 
values  of  V  and  K,  it  appears  that  the  analytical  model  can 
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Figure    (4)    -  Formaldehyde  Concentrations    in  an   Idealized 
Two-dimensional   Columnar   Photo-reactor 
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be  significantly  improved  (see  Figure (4))  by  employing  a 
higher  intensity  mercury  lamp  (increasing  K)  and/or  by 
decreasing  the  flow  velocity. 

Increasing  the  intensity  of  the  lamp  is  a  viable 
option.   There  are  currently  bulbs  available  which  are 
capable  of  over  1000  times  the  power  of  the  mercury  vapor 
lamps  which  are  used  in  the  photolytic  reactor  of  Figure 
(2)  . 

Since  a  decrease  in  velocity  beyond  some  practical 
value  will  only  lead  to  a  decrease  in  the  amount  air 
treated,  increasing  the  bulb  intensity  is  clearly  the  best 
practical  alternative. 
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Chapter  7 
DEVELOPMENT  OF  THE  NUMERICAL  MODEL 

The  numerical  photolytic  reactor  model  under 
consideration  predicts  the  concentrations  in  a  rectangular 
area  around  a  typical  bulb  in  the  reactor  of  Figure  (2) . 
Since  very  little  reduction  in  formaldehyde  is  expected  in  a 
small  21.2  x  6.6  cm  domain  utilizing  the  actual  rate 
constant  calculated  in  chapter  5.1.1,  the  rate  constant  has 
been  increased  10,000  times  to  an  attainable  0.35  (1/s)  to 
achieve  some  discernable  results  within  two  decimal  places. 
Although  this  bulb  intensity/reaction  rate  increase  is 
presently  impractical  from  a  cost  standpoint  (very  expensive 
bulbs  are  required)  it  greatly  improves  the  analysis  and 
interpretation  of  the  numerical  results  for  the  small  flow 
domain  under  consideration.   The  following  data  describe  the 
parameters  of  the  model: 

L,  length  of  the  domain  in  x  direction  =  .2121  m 
M,  number  of  nodes  in  the  x  direction  =  100 
h,  size  of  increments  in  the  x  direction  =  0.002143  m 
H,  length  of  the  domain  in  the  y  direction  =  .0664  m 
N,  number  of  nodes  in  the  y  direction  =  31 
h,  size  of  increments  in  the  y  direction  =  0.002143  m 
Uo,  inlet  velocity  =  2.0  m/s  (const,  across  entire 
area) 
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Co,  inlet  formaldehyde  concentration  =50.0  PPM  (  "  ) 
Bulb  diameter  =  .03  m 
K,  overall  rate  constant  at  bulb  =  .35  1/s 

D,  diffusion  coefficient  =  0.00002  m2/s 

£ ,  solution  convergence  criteria  =  0.000001  PPM 

(maximum  difference  between  successive  concentration 
calculations  at  each  point  at  a  given  time) 
At,  time  increment  =  0.01  s 


The  computer  model  uses  a  finite  difference  technique 
in  conjunction  with  a  Gauss-Seidel  iterative  method  to 
calculate  the  concentration  field  throughout  the  reactor  at 
selected  time  intervals.   The  mass  conservation  equation  is 
restated  as  follows: 

-voVC  +  DV2C  +  Q  =  dC/dt  [31] 

Since  the  convection  terms  dominate  Equation  [31],  an 
upwind  difference  method  has  been  incorporated  into  the 
above  equation  to  improve  the  accuracy  of  the  results. 
Simply  stated,  the  upwind  technique  allows  the  concentration 
at  any  point  i, j  to  be  calculated  based  on  points  upwind 
since  concentration  changes  are  much  more  closely  related  to 
the  fluid  flow  than  they  are  to  diffusion  and  photolytic 
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effects.   Hence,  differencing  in  the  upwind  direction  is 
utilized  for  the  convection  terms.   For  the  sake  of 
simplicity,  central  differencing  is  used  for  the  diffusion 
terms.   The  convection  term  of  Equation  [31]  will  now  be 
developed  in  terms  of  finite  differences. 

The  equation  which  follows,  Equation  [32],  is  an  upwind 
difference  expression  which  is  valid  when  both  the  fluid 
velocity  components,  vx  and  v  are  positive. 

-VoVC  =  -(vx/h)  (Cifj-Ci_1/j)-(vy/h)  (C±rrC±f1j_1)     [32] 

When  both  the  fluid  velocity  components  are  negative: 

-VoVC  =  -(vx/h)  (C±+lf  j-Cif  j)-(vy/h)  (Ci/j+1-Ci/:j)  [33] 

An  expression  derived  from  Equations  [32]  and  [33]  above 
which  will  always  yield  an  upwind  difference  for  any 
combination  of  velocity  component  directions  can  be  written 
as  follows: 

-voVC  =  -(l/2h)(  (|vx|+  vx)  (C±/  j-Ci_1/  j)  + 

(|vy|+  vy)  (Cifj-Cifj_1)  +  (vx-  |vx|)(C1+lfj-Clfj)  + 

(vy-  lvyD  (ci,j  +  l"ci,j)'  t3^ 

While  the  above  equation  may  appear  unnecessarily  lengthy, 
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it  greatly  increases  the  speed  of  numerical  calculations  as 
a  decision  does  not  have  to  be  made  (FORTRAN  IF  statement) 
at  each  point  in  order  to  apply  the  correct  equations  based 
on  the  velocity  component  directions.   Without  condition 
statements  in  the  calculation  of  the  convective  terms,  much 
of  the  program  can  be  vectorized  to  greatly  reduce  the 
computation  time. 

In  terms  of  central  differences,  the  diffusion  term  can 
be  written: 


DV2C  =   (D/h*>fc1+1  ,  ♦  Cl.lfj  - 
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4ci,j  +  ci,j+l  +  ci,j-l'  [35] 


The  remaining  terms  are 


Q.      =   -KC4  i  [36] 

dC/dt  =  (1/At)  (Cj  ,n  -  C,  in~1)  [37] 


The  superscripts  n  and  n-1  above  refer  to  the  present 
value  and  the  value  at  the  previous  time,  respectively.   It 
is  important  to  note  that  a  backward  difference  scheme  has 
been  employed  for  the  time  dependant  solution  of  the 
governing  equation.   Although  all  terms  on  the  left  side  of 
the  governing  equation  (Eqn.  [31])  have  been  developed  using 
concentration  values  at  the  present  time,  the  superscript  n 
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has  been  omitted  for  clarity.   The  backward  difference 
method  results  in  an  implicit  set  of  equations  which  can  be 
solved  by  the  Gauss-Seidel  iterative  method  to  find  the 
concentration  at  each  point  in  the  reactor  at  a  particular 
point  (iteration)  in  time,  n.   This  method  may  seem 
cumbersome  in  comparison  to  the  forward  difference  method 
which  can  explicitly  calculate  the  concentrations  at  a 
future  point  in  time  without  iteration.   However,  the 
backward  difference  method  has  higher  stability  which  allows 
a  much  larger  time  step  to  be  used. 

Substituting  Equations  [33]  through  [37]  into  Equation 
[28]  and  solving  for  C^  a    yields  the  following  iterative 
equation: 


(l/2)f  (|vj+  vv) 


I  lvxl+  lvyl+  (4D/h)  +  hK  +  (h/At) 

Ci-i,jk+1  +  (|vy|+  VCi.J-lk+1  +  (|vxl"  vx)Ci+l,jk  + 

<  I  vy  I  "  vy)Ci,j  +  lkl  +  <D/h>fci  +  l,jk  +  Ci-l,jk+1  + 
ci,j+lk  +  ci,j-lk+1'  +  (h/At)Ci^jn-1)  [38] 

The  superscripts  k  and  k+1  in  equation  [38]  represent 
concentration  values  calculated  at  the  past  and  present 
iterations,  respectively.   Note  that  the  Ci_1  j  and  C^  j_x 
terms  have  the  superscript  k+1  whereas  the  Ci+1  a   and 
Cj_  j+1  terms  have  the  superscript  k.   This  is  due  the 
iterative  "search  pattern"  which  sweeps  from  j=l  to 
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N  (y=0  to  y=H)  for  each  i  until  i=M  (x=L)  is  reached.   This 
particular  sweep  pattern  in  conjunction  with  the  Gauss- 
Seidel  method  makes  use  of  the  latest  values  (k+1)  for 
Ci_1  j  and  C-l  j_x  in  calculating  the  updated  Ci  jk+1  value. 
The  boundary  conditions  of  the  flow  domain  are  set  as 
follows : 

A.  Inlet  -  C±  =  const.  =  Co  (i=l; l<j>25)  [39] 

B.  Top  and  bottom  surfaces  -  dC/dy  =  0  (symmetry)  or, 

ci,j+l  =  ci,j-l  (l*i*100;j-l,25)  [40] 

C.  At  bulb  -  dC/dx,  dC/dy  =  0  (solid  boundary)  or, 

Ci+1  j  =  Ci_1  .j  (see  program,  Appendix  A)      [41] 
ci  j+1  =  ci  j-1  (see  program,  Appendix  A)      [42] 

D.  Exit  -  dC/dx  =  0  (minimal  change  at  exit)  or, 

ci,j  ■  ci-l,j  (i=100;l<j>25)  [43] 

The  computer  program  included  in  Appendix  A  uses 

Equation  [38]  (including  the  relationships  of  Equations  [39] 

through  [43]  when  required)  to  calculate  and  update  new 

concentration  values  at  all  grid  points  in  the  domain  per 

each  iteration.   The  iterations  converge  to  the  solution  for 

a  given  time  step  when  the  present  and  previous 

concentration  values  for  each  point  in  the  entire  domain 

k+1  _c   k| 


differ  by  less  than  an  epsilon  value  (i.e., 
<  £  for  all  points  in  the  domain) . 


r. 


Figure  (5)  is  a  graphical  representation  of  the  program 
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Figure  (5)  -  Formaldehyde  Concentrations  Within  the 
Numerical  Photo-reactor  Model 
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output  for  selected  values  of  time.   The  solution  at  t  =  10 
sec  is  considered  to  be  a  steady  state  solution.   This  was 
determined  by  comparing  the  solutions  at  t  =  10  sec  and  t  = 
100  sec  which  did  not  differ  to  two  decimal  places  for  the 
vast  majority  of  points  in  the  domain.   In  addition,  to 
ensure  that  a  time-accurate,  steady  state  solution  was 
obtained,  the  computer  program  was  run  at  At=0.005  s  (half 
the  original  time  step)  out  to  a  10  sec  solution.   The 
results  of  this  run  did  not  differ  from  the  original  10  sec 
solution  to  over  4  decimal  places  for  the  entire  domain. 

The  numerical  concentration  predictions  for  the  "start- 
up" of  the  reactor  make  sense  from  a  practical  standpoint. 
The  smallest  concentrations  are  seen  at  the  bulb  surface  and 
wake  regions  due  to  the  low  fluid  velocities  there.   In 
addition,   these  regions  of  lesser  concentration  migrate 
further  downstream  with  increased  time  as  they  are  swept 
into  the  bulb  wake  region. 

The  steady  state  isopleths  also  appear  reasonable  in 
form.   As  expected,  these  concentration  contours  generally 
form  the  same  shape  as  the  streamlines  depicted  in  the 
FLUENT  output  due  to  the  dominance  of  the  convection  terms. 
However,  the  isopleth  shapes  differ  from  the  streamlines  due 
to  diffusive  transport  and  photolysis.   In  the  next  chapter 
it  will  be  shown  by  comparison  with  the  simplified  flowfield 
model  that  the  magnitude  of  formaldehyde  reductions  shown  in 
Figure  (5)  are  reasonable  predictions. 
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Chapter  8 
ANALYSIS  OF  RESULTS 

In  order  to  obtain  a  general  performance  check  on  the 
numerical  solution,  the  simplified  flow  model  developed  in 
Chapter  5  can  be  employed  by  inputting  the  numerical  model 
parameters  discussed  in  Chapter  7  into  Equation  [30]  as 
follows: 

Co  =  50  (PPM) 
K  =  0.35  (1/s) 
v  =  2.0  (m/s) 
x  =  0.2121  (m) 

Substituting  values  into  equation  [30]  yields: 

C  =  Co  exp{-Kx/v} 

=  50.0PPM  exp{-0. 35s"1 (0.2121m) /2.0m/s} 
=  48.17  PPM 

The  above  result  compares  favorably  with  the  4  9.37  PPM 
average  steady  state  concentration  leaving  the  numerical 
reactor  model  which  was  found  through  a  mass  average  taken 
across  the  exit  as  follows: 

Cexit,avg.  =  ?>V±C±/I,v±  [45] 
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In  the  next  paragraph  the  analytical  and  numerical  results 
are  compared  in  terms  of  destruction  efficiencies. 

An  additional  comparison  between  the  numerical  and 
simplified  flow  field  models  is  obtained  by  introducing 
parameters  which  define  the  destructive  efficiency  of  the 
photolytic  reactor.   Since  it  is  the  objective  of  the 
reactor  model  to  eliminate  formaldehyde,  the  efficiency  of 
the  entire  reactor  section  can  be  defined  as  follows: 

E  =  x  "  (cexit,avg./cinlet>  [46] 

In  addition  to  the  above  equation  it  is  also  useful  to 
define  an  efficiency  which  is  based  on  the  bulb  height 
alone.   This  term  is  defined  as  follows: 


EB  -  1  -  (Cexit/avg-/Cinlet)  [47] 


Equation  [47]  calculates  an  efficiency  based  on  a  collimated 
region  with  a  height  equal  to  the  radius  of  the  bulb.   In 
Equation  [47]  the  term  Cexit,avg. /cinlet  is  taken  from  the 
centerline  of  the  reactor  up  to  the  bulb  radius.   Table-1 
compares  the  efficiency  of  the  steady  state  numerical 
simulation  of  the  photolytic  reactor  section  with  that  of 
the  simplified  flow  field  model. 
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Numerical  Model 
(PPM) 

Simplified  flow 

Field  Model 

(PPM) 

C  (average  across 
entire  reactor 
section) 

49.37 

48.17 

C  (average  up  to 
bulb  radius) 

47.37 

48.17 

E  (based  on  entire 
reactor  cross 
section) 

1.26% 

3.64% 

EB  (based  on  bulb 
width) 

5.26% 

3.64% 

Table-1   Comparison  of  Numerical  and  Simplified  Flow  Models 


For  this  set  of  conditions  it  appears  that  the 
simplified  flow  field  model  is  more  efficient  than  the 
numerical  reactor  model  when  the  efficiency  is  based  on  the 
entire  reactor  cross  section.   However,  in  examining  the 
efficiencies  based  on  the  bulb  cross  section,  the  numerical 
reactor  model  which  incorporates  the  actual  flow  conditions 
is  superior  to  the  simplified  flow  field  model  in  its 
ability  to  destroy  the  formaldehyde.   This  result  leads  one 
to  believe  that  a  reactor  containing  a  combination  of  bulbs 
which  allows  the  incoming  flow  to  "see"  bulbs  across  the 
entire  reactor  area  will  out-perform  the  idealized  columnar 
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reactor. 

The  advantage  of  using  multiple  bulbs  is  best 
illustrated  by  considering  a  single  bulb  reactor  of  infinite 
height.   In  such  a  reactor  the  UV  bulb  would  have  little 
effect  on  the  VOC  concentrations  in  the  far  stream.   As 
shown  in  Figure  (5) ,  the  bulb  has  the  most  influence 
(concentration  reductions)  in  a  columnar  region  which  is  one 
to  two  bulb  diameters  in  height.   In  other  words,  the 
incoming  flow  which  exists  outside  this  columnar  region 
tends  to  escape  untreated.   Therefore,  an  arrangement  of 
bulbs  which  guarantees  that  all  the  incoming  flow  will 
encounter  a  bulb  cross  section  (ie.,  a  staggered  stack  of 
bulbs/columnar  regions  in  parallel  with  no  "escape"  regions 
in  between)  would  have  a  destructive  efficiency  of  no  less 
than  that  given  by  Equation  [47] . 

In  addition  to  the  general  advantages  of  a  parallel 
bulb  arrangement  as  described  above,  large-scale  destructive 
efficiency  improvements  for  both  series  and  parallel 
arrangements  are  predicted  by  drawing  an  analogy  between 
particle  collection  and  VOC  destruction  in  a  photolytic 
reactor.   For  such  an  analogy  the  multi-bulb  efficiency  is 
given  by  the  following  (Heinsohn,  1991) : 


'overall 


=  1  -  (1  -  ER)N  [48] 


where,  N  is  the  number  of  bulbs 
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A  twenty-five  (N  =  25)  bulb  arrangement  with  single  bulb 
efficiencies  as  given  by  Equation  [47]  (EB  =  5.26%  in  this 
case)  would  have  an  overall  destructive  efficiency  as 
follows : 


Overall  =  1  "  d  "  EB)N   =  1  -  (1  -  0.0526) 25 
=  74.1% 


Although  the  particle  collector/photolytic  reactor 
analogy  drawn  here  predicts  very  optimistic  results,  it  will 
be  necessary  to  confirm  this  relationship  with  a  multi-bulb 
numerical  model  and/or  experimental  testing  to  achieve  a 
high  degree  of  confidence. 
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Chapter  9 
CONCLUSIONS 

This  analysis  clearly  demonstrates  the  influence  of 
actual  flow  conditions  on  the  effectiveness  of  a  photolytic 
reactor.   Despite  the  simplifying  (but  reasonable) 
assumptions  made  in  development  of  the  numerical  model,  it 
provides  valuable  insight  into  the  inner  workings  of  the 
actual  photolytic  reactor.   The  graphical  results  of  Figure 
(5)  show  that  the  velocity  gradients,  recirculation  and 
mixing  which  occur  within  the  reactor  contribute  greatly  to 
the  overall  reduction  of  reactants.   Since  recirculation  and 
mixing  greatly  increase  the  reduction  of  formaldehyde  in  the 
reactor,  one  would  expect  that  unsteadiness  in  the  flow 
field  would  tend  to  improve  the  results  (noting  that 
unsteadiness  exists  in  the  actual  flow  field  under 
consideration  due  to  vortex  shedding) .   In  addition,  one 
would  expect  that  an  optimum  Reynolds  number,  Re  could  be 
found  (input  velocity)  which  would  maximize  the  benefits  of 
mixing  while  still  providing  an  acceptable  level  of  air 
treatment.   Therefore,  future  numerical  optimization  studies 
should  include  both  an  unsteady  flow  field  analysis  and  a 
comparison  of  results  for  different  Reynolds  numbers. 

Although  the  simplified  flow  field  model  developed  in 
this  report  can  be  used  to  obtain  a  rough  estimate  of 
reactor  performance,  it  is  unable  to  predict  the  effects  of 
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the  actual  flow  field  which  were  found  to  be  of  great 
importance  through  an  analysis  of  the  numerical  model 
results.   This  fact  reinforces  the  need  to  use  a  numerical 
model  in  the  optimization  of  a  photolytic  reactor  design. 

Comparison  of  efficiencies  between  the  simplified  flow 
field  and  numerical  models  provides  some  highly  instructive 
results.   Surprisingly,  though  the  light  intensity  in  the 
numerical  model  is  more  than  five  times  less  than  that  of 
the  simplified  flow  field  model  for  a  large  portion  of  the 
domain  (due  to  the  radial  light  attenuation  from  the  lamp 
surface) ,  the  numerical  model  provides  a  comparable 
formaldehyde  reduction  performance  when  the  efficiency  is 
based  on  the  entire  reactor  cross  section.   In  addition,  it 
appears  that  increasing  the  number  of  bulbs  will  greatly 
improve  overall  reactor  effectiveness  as  indicated  by  a 
comparison  of  efficiencies  based  on  the  bulb  width.   Due  to 
the  optimistic  results  obtained  in  comparing  the  numerical 
and  simplified  flow  field  models,  the  next  logical  step  for 
future  research  in  this  area  should  include  the  modeling  of 
larger  flow  domains  which  contain  a  larger  quantity  of 
bulbs. 

This  undertaking  also  showed  that  additional  factors 
can  readily  be  taken  into  consideration  in  future  reactor 
modeling.   The  computer  program  of  Appendix  A  and  techniques 
described  in  this  analysis  can  be  extended  to  include 
additional  chemical  reactions  and  other  factors  such  as 
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Beer-Lambert  light  attenuation  by  reactants  which  will 
result  in  a  numerical  model  of  much  greater  accuracy. 
Although  this  work  is  of  limited  scope  for  practical 
reasons,  expansion  of  this  analysis  including  the  use  of  a 
high-speed  supercomputer  could  achieve  an  accurate  numerical 
simulation  of  the  entire  photolytic  reactor  of  Figure  (2) . 

Despite  the  good  correlation  between  the  numerical  and 
simplified  flow  field  models  developed  in  this  report,  a 
high  degree  of  confidence  can  only  be  achieved  through 
laboratory  experiments  whose  results  can  be  compared  to  the 
numerical  ones.   Therefore,  future  work  should  include 
experimental  verification.   This  will  be  especially 
important  when  additional  reactants  and  other  factors  are 
introduced.   In  the  future,  the  benefits  of  a  closed  form 
solution  as  experienced  in  this  development  (simplified  flow 
field  model)  will  not  be  possible  when  modeling  a  highly 
complex  reactor. 
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Appendix  A 
COMPUTER  PROGRAM 


50 

C  SET  UP  VARIABLES 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

DIMENSION  VX(100,33),VY(100,33) ,FMC(100,33) ,FMCOLD (100, 33) 

X0=4.286 

Y0=0.0 

R0=1.50 

M=100 

N=31 

H=0. 002143 

EPSI=0. 0000001 

DT=0.01 

ITER=1 

ITERT=1 

ITEST=1 

ITERMAX=200 

ITERTMAX=10 
C  ENTER  FORMALDEHYDE  CONSTANTS 

FMRC0=0.35 

FMD=0. 00002 

FMC0=50.0 
C  READ  IN  FLUENT  GENERATED  VELOCITIES 

OPEN (UNIT=5,FILE=' SINGLEBULBLIS.DAT' , STATUS=' OLD' ) 

K=l 

KK=1 

KKK=10 

2  READ(5,5)   ( (VX (I, 33- J) , I=KK, KKK) , J=l ,N) 
K=K+1 

KKK=K*10 

KK=(K*10) -9 

IF(K.LE.IO)  GO  TO  2 

K=l 

KK=1 

KKK=10 

3  READ(5,5)   ( (VY (I, 33-J) , I=KK, KKK) , J«1,N) 
K=K+1 

KKK=K*10 
KK=(K*10) -9 
IF (K.LE.10)  GO  TO  3 
CLOSE (UNIT=5) 
5  FORMAT (5X,10E12. 4) 
C  SET  THE  INITIAL  FORMALDEHYDE  CONCENTRATIONS 
DO  10  1=1, M 
DO  10  J=l,N+2 
FMC(I,  J)=FMC0 
FMCOLD(I,  J)=FMC0 
10  CONTINUE 
C  CALCULATE  THE  CONCENTRATIONS  AT  EACH  POINT  USING  GAUSS-SEIDEL  METHOD 
15  DO  50  1=2, M-l 
DO  50  J=2,N+1 
TEMP=FMC(I,  J) 
C   CALCULATE  THE  PHOTOLYTIC  RATE  CONSTANT 

FMRC=FMRC0*R0/ ( (100. *H* (1-1) -X0) **2 . + (100 . *H* (J-2) ) **2.+ 
$0.00000000001) **0.5 
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C  APPLY  BOUNDARY  CONDITIONS  (BC'S)  AT  THE  BULB  (J=2, 1=14 ; J=4 , 1=15 ; 
C  J=6,I=16, J=7, 1=17, J=8, 1=19) 

IF( ( (J.EQ.2)  .AND.  (I.EQ.14) )  .OR.  ( (J.EQ.4)  .AND.  (I.EQ.15) )  .OR. 
$ ( (J.EQ.6)  .AND.  (I.EQ.16) )  .OR.  ( (J.EQ.7)  .AND.  (I.EQ.17) )  .OR. 
$ ( (J.EQ.8) .AND. (I.EQ.19) ) ) 

$FMC(I, J)=(l/ (4.*FMD/H**2.+FMRC+1./DT) ) * 

$ ( (FMD/H**2.) *(2.*FMC(I-1, J)+2.*FMC(I, J+l) ) +FMCOLD (I, J) /DT) 
C  APPLY  BOUNDARY  CONDITIONS  (BC'S)  AT  THE  BULB  ( J=2 , 1=28 ; J=4 , 1=27 ; 
C  J=6, 1=26, J=7, 1=25, J=8, 1=23) 

IF ( ( (J.EQ.2)  .AND.  (I.EQ.28) )  .OR.  ( (J.EQ.4)  .AND.  (I.EQ.27) )  .OR. 
$ ( (J.EQ.6)  .AND.  (I.EQ.26) )  .OR.  ( (J.EQ.7)  .AND.  (I.EQ.25) )  .OR. 
$ ( (J.EQ.8) .AND. (I.EQ.23) ) ) 

$FMC(I, J)=(l/ (4.*FMD/H**2.+FMRC+1./DT) ) * 

$ ( (FMD/H**2.) * (2.*FMC(I+1, J)+2.*FMC(I, J+l) ) +FMCOLD (I, J) /DT) 
C   SKIP  AREAS  INSIDE  THE  BULB 

IF( ( ( (J.GE.2) .AND. (J.LE.3) ) .AND. ( (I.GE.15) .AND. (I.LE.27) ) ) .OR. 
$ ( ( (J.GE.4) .AND. (J.LE.5) ) .AND. ( (I.GE.16) .AND. (I.LE.2  6) ) ) .OR. 
$ ( (J.EQ.6) .AND. ( (I.GE.17) .AND. (I.LE.25) ) ) .OR. ( (J.EQ.7) .AND. 
$ ( (I.GE.18) .AND. (I.LE.24) ) ) .OR. ( (J.EQ.8) .AND. ( (I.GE.20) .AND. 
$ (I.LE.22) ) ) )  GO  TO  50 
C   CALCULATION  FOR  GENERAL  POINTS 
IF( (J.GE.9) .OR. 
$ ( ( (I.LE.13) .OR. (I.GE.2  9) ) .AND. (J.EQ.2) ) .OR. 

$ ( ( (I.LE.14) .OR. (I.GE.28) ) .AND. ( (J.GE.3) .AND. (J.LE.4) ) ) .OR. 
$ ( ( (I.LE.15) .OR. (I.GE.27) ) .AND. ( (J.GE.5) .AND. (J.LE.6) ) ) .OR. 
$ ( ( (I.LE.16) .OR. (I.GE.2  6) ) .AND. (J.EQ.7) ) .OR. 
$ ( ( (I.LE.18) .OR. (I.GE.24) ) .AND. (J.EQ.8) ) ) 

$FMC(I, J)=(l/ ( (ABS(Vx(I, J) )+ABS (VY(I, J) ) +4 . *FMD/H) /H+FMRC+1 . /DT) ) * 
$  (  (  (ABS(VX(I,  J)  )+VX(I,  J)  )  *FMC(I-1,  J)  +  (ABS  (VX  (I ,  J)  )  -VX(I,  J)  )  * 
$FMC(I  +  1,  J)  +  (ABS(VY(I,  J)  )+VY(I,  J)  )  *FMC(I,  J-l)  + 
$  (ABS(VY(I,  J)  )  -VY(I,  J)  )  *FMC(I,  J+l)  )  /  (2 . *H) + (FMD/H**2 . )  * 
$  (FMC(I-1,  J)+FMC(I  +  1,  J)+FMC(I,  J-1)+FMC(I,  J+l)  ) +FMCOLD  ( I ,  J)  /DT) 
IF (ABS (FMC(I, J) -TEMP) .GT.EPSI)  ITEST=ITEST*0 . 
50  CONTINUE 
C  ASSIGN  VALUES  REQUIRED  DUE  TO  BOUNDARY  CONDITIONS 
DO  52  1=2, M-l 
FMC(I,1)=FMC(I,3) 
FMC(I,33)=FMC(I,31) 
52  CONTINUE 

FMC(20,8)=FMC(20,10) 
FMC(21,8)=FMC(21,10) 
FMC(22,8)=FMC(22,10) 
FMC(18,7)=FMC(18,9) 
FMC(24,7)=FMC(24,9) 
FMC(16,5)=FMC(14,5) 
FMC(26,5)=FMC(28,5) 
FMC(15,3)=FMC(13,3) 
FMC(27,3)=FMC(29,3) 
C   APPLY  BC'S  AT  END  (DC/DX=0) 
DO  55  J=l,N+2 
55  FMC(M,  J)=FMC(M-1,  J) 

IF(ITER.GE.ITERMAX)  GO  TO  60 
IF(ITEST.EQ.l)  GO  TO  60 
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ITER=ITER+1 

ITEST=1 

GO  TO  15 
60  IF(ITERT.EQ.ITERTMAX)  GO  TO  64 

DO  62  1=2, M 

DO  62  J=2,N+1 
62  FMCOLD(I,  J)=FMC(I,  J) 

ITER=1 

ITERT=ITERT+1 

ITEST=1 

GO  TO  15 
WRITE  OUT  RESULTS 

64  OPEN (UNIT=6,FILE=' NEWS INGLE. OUT' , STATUS='NEW' ) 
OPEN  (UNIT=8,FILE='  SINGLEPLOT.DAT'  , STATUS=' NEW'  ) 
DO  100  1=15,27 

DO  100  J=2,8 
ASSIGN  VALUES  TO  ACHIEVE  A  DECENT  CONTOUR  PLOT  USING  GNUPLOT 
100  IF( ( ( (J.GE.2) .AND. (J.LE.3) ) .AND. ( (I.GE.15) .AND. (I.LE.27) ) ) .OR. 
$ ( ( (J.GE.4) .AND. (J.LE.5) ) .AND. ( (I.GE.16) .AND. (I.LE.26) ) ) .OR. 
$ ( (J.EQ.6)  .AND.  ( (I.GE.17)  .AND.  (I.LE.25) ) )  .OR.  ( (J.EQ.7)  .AND. 
$ ( (I.GE.18) .AND. (I.LE.24) ) ) .OR. ( (J.EQ.8) .AND. ( (I.GE.20) .AND. 
$ (I.LE.22) ) ) )  FMC(I,J)=30. 
SUM=0. 
CALCULATE  EXIT  AVERAGES 
SUMD=0. 
SUM1=0. 
SUM1D=0 . 
DO  80  J=2,N+1 
SUMD=SUMD+VX ( 1 0  0 ,  J) 

80  SUM=SUM+FMC(100, J) *VX(100,  J) 
SUM=SUM/SUMD 

DO  81  J=2,9 
SUM1D=SUM1D+VX ( 1 0 0 ,  J) 

81  SUM1=SUM1+FMC(100, J) *VX(100,  J) 
SUM1=SUM1/SUM1D 

K=l 

KK=1 

KKK=10 

WRITE (6,71)  ITER, ITERT , SUM, SUM1 , DT , EPS I 

65  WRITE (6, 70)   ( (FMC (I, 33-J)  ,  I=KK,KKK) , J=l, N) 
K=K+1 

KKK=K*10 
KK=(K*10)  -9 
IF(K.LE.IO)  GO  TO  65 
DO  67,  J=2,N+1 
DO  66,  1=1, M 

66  WRITE  (8, 68)  FMC (I, J) 
WRITE (8, 69) 

67  CONTINUE 

68  FORMAT  (F12. 8) 

69  FORMAT (IX) 

70  FORMAT (10F12. 8) 

71  FORMAT (216, 4E12 . 4) 
STOP 

END 
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Appendix  B 
ATMOSPHERIC  CHEMISTRY 
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WELL-MIXED  MODEL  OF 
PHOTO-OXIDATION  AND  RADICAL  OXIDATION  PROCESSES 

Contents: 

Photolysis  Rates  for  337  nm (laser) 2 

Photolysis  Rates  for  253  &  185  nm. (UV  bulb) 3 

Photolysis  Rates  for  laser  &  UV  bulb (combo) 4 

Radical  Chemistry 5 


Reactions: 

Formaldehyde. 6 

Methanol 
Carbon  Dioxide 

Water  Vapor. 7 

Ozone 


Atomic  Oxygen 8 

NO 9 

H2O2 

N02 10 

HONO 

Hydroxyl  Radicals..(OH) 1 1 

Hydroperoxyl  Radicals..(H02) 1 2 

Carbon  Monoxide  (CO) 

Atomic  Hydrogen.. (H) 13 

Formyl  Radical  (HCO) 

Atmospheric  Gases 14 

References. 15 

Note: 

R  represents  a  reaction;  P  represents  a  photolytic  reaction;  KR 

represents  a  rate  constant;  KP  represents  a  photolytic  rate  constant. 
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Photolysis: 

X  for  the  laser  experiment  will  be  337  nm 
input  file  "RATE" 
Ephoton=590*10-21  J/Photon 

Actinic  Flux=20*1015  Photons/(cm2*s) 


HCHO+hv-*H*+HCO* 
HCH0+hv-COH2 
CH30H+hv  —  negligible 

Os+hv  —  O^D)-^ 
03+hv  —  0(3p)+02 
02+hv-*0(1D)+0(3P) 
02+hv-0(3p)+0(3p) 

Q2+hv-0(lD)+0(lD) 

H20z+hv-*0H*+0H* 
H202+hv-H204O(1D) 
H202+hv  -»H*+H02* 
H202+hv^2H*+02 

N02+hv-N0+0(lD) 
N02+hv-*NO+0(3p) 

H20+hv-H2+0(3p) 
H20+hv— H*+OH* 
H20+hv—H2+0(1D) 

C02+hv-*CO0(lD) 
CCb+hv  — C04O(3p) 

HONOhv-NOOH* 

HONOhv— H*-rN02 
HONO+hv  -*HN0*+0(3P) 


\<300nm 

4.08E-5  s'1 

[P1A] 

Atkinson 

>^300nm 

2.97E-4  s"1 

[P1B] 

Atkinson 

0 

[P2] 

>^41  lnm 

0 

[P3A] 

Atkinson 

X<1180nm 

0,  (O  =  0) 

[P3B] 

Atkinson 

>.<175nm 

0 

[P4A] 

Atkinson 

X<242nm 

0 

[P4B] 

Johnson/ 
Atkinson 

X<137nm 

0 

[P4C] 

Johnson/ 
Atkinson 

>.<557nm 

1.40E-5  s-1 

[P5A] 

Atkinson 

X<359nm 

0,(0  =  0) 

[P5B] 

Atkinson 

>.<324nm 

0,(4>=0) 

[P5C] 

Atkinson 

X<213nm 

0 

[PSD] 

Atkinson 

Xs244nm 

0 

[PGA] 

Atkinson 

>.<398nm 

7.08E-3  s"1 

[P6B] 

Atkinson  * 

X<243nm 

0 

[P7AJ 

Atkinson 

X<239nm 

0 

P7B] 

Atkinson 

X<176nm 

0 

negligible 
negligible 

[P7C] 

Atkinson 

X<591nm 

3.23E-3S-1 

[P50A] 

Atkinson 

X<367nm 

0,  (0  =  0) 

[P50B] 

Atkinson 

X<283nm 

0 

[P50C] 

Atkinson 

Photolysis: 
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X  for  the  UV  Bulb  experiment  will  be  253  nm  and  185  nm 
input  file  "BULB" 
Ephoton=5  90*  10-21  j/Photon 

Actinic  Flux:  J(253)=1.14*1016  Photons/(cm2*s) 
J(185)=3.52*10l4  Photons/(cm2*s) 


HCHOhv-*H*+HC0* 

X<300nm 

1.30E-5S"1 

[P1A] 

Atkinson 

HCHO+hv-COH2 

X*300nm 

2.16E-5  s"1 

[P1B] 

Atkinson 

CH30H+hv  -*  negligible 

[P2] 

Q3+1*  — O^DJ+Qz 

>cs411nm 

1.12E-1  s'1 

[P3A] 

Atkinson 

03+hv^O(3P)+02 

>.<1180nm 

1.32E-2  s'1 

-  [P3B] 

Atkinson 

02+hv—0(1D)+0(3p) 

X<175nm 

0 

[P4A] 

Atkinson 

02+hv-0(3p)+0(3p) 

X<242nm 

1.23E-9  s'1 

[P4B] 

Johnson/ 
Atkinson 

07+hv->0(iD)+0(iD) 

X<137nm 

0 

[F4CJ 

Johnson/ 
Atkinson 

H202+hv-*OH*+OH* 

X<557nm 

8.38E-4  s"1 

[P5A] 

Atkinson 

H202+hv  — H20+0(1D) 

X<359nm 

0,(0  =  0) 

[P5B] 

Atkinson 

H202+hv  -*H*+H02* 

X<324nm 

0,  (O  =  0) 

[P5C] 

Atkinson 

H202+hv-*2H*+02 

X<213nm 

0,(0  =  0) 

[PSD] 

Atkinson 

N02+hv-NOO(1D) 

Xs244nm 

0,(0  =  0) 

[P6A] 

Atkinson 

N02+hv  — N0+0(3P) 

X<398nm 

2.47E-4  s'1 

[P6B] 

Atkinson 

H20+hv-*H2+0(3p) 

X<243nm 

0,  (0  =  0) 

[P7A] 

Atkinson 

H20+hv^H*+0H* 

X<239nm 

1.09E-5  s'1 

[P7B] 

Atkinson 

H20hv-*H2+0(1D) 

X<176nm 

0 

[P7C] 

Atkinson 

C02+hv-CO0(lD) 

negligible 

C02+hv— COO(3P) 

negligible 

HONO+hv-NO+OH* 

X<591nm 

1.69E-3  s"1 

[P50A] 

Atkinson 

HONO+hv—  H*+N02 

X<367nm 

0,  (0  =  0) 

[P50B] 

Atkinson 

HONO+hv  — HNO*+0(3P) 

X<283nm 

0,(0  =  0) 

[P50C] 

Atkinson 
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Photolysis: 


>w  for  the  UV  Bulb  and  the  laser  experiment  will  be  337,  253,  and  185  nm 
input  file  "COMBO" 
Ephoton=590*10-21  J/Photon 


Actinic  Flux: 


J(337)=20.0*1015  Photons/(cm2*s) 
J(253)=l. 14*1016  Photons/(cm2*s) 
J(185)=3.52*1014  Photons/(cm2*s) 


HCHCM-hv-H*+HC0* 

X<300nm 

1.30E-5s  -1 

[P1A] 

Atkinson 

HCHOhv-COfHz 

Xa300nm 

3.86E-5  s"1 

[P1B] 

Atkinson 

CH3OH+I1V  -+  negligible 

[P2] 

03+hv— O^D)-^ 

X^411nm 

1.12E-1  s"1 

[P3A] 

Atkinson 

03+hv  —  0(3P)+02 

X<1180nm 

1.32E-2  s'1 

[P3BJ 

Atkinson 

02+hv-*0(lD)+0(3p) 

>.<175nm 

0 

[P4A] 

Atkinson 

02+hv-0(3p)+0(3p) 

X<242nm 

1.23E-9  s"1 

[P4BJ 

Johnson/ 
A  rki  n  son 

02+hv-0(1D)+0(lD) 

>w<137nm 

0 

[P4C] 

Johnson/ 
Atkinson 

H202+hv-*OH*+OH* 

X<557nm 

8.52E-4  s'1 

[P5A] 

Atkinson 

H202+hv-H204O(lD) 

>*<359nm 

0,  (0  =  0) 

[P5B] 

Atkinson 

H202+hv—H*+H02* 

>v<324nm 

0,  (O  =  0) 

[P5C] 

Atkinson 

H202+hv^2H*+02   " 

>*<213nm 

0,  (0  =  0) 

[P5D] 

Atkinson 

NQ2+hv  — N0+0(1D) 

Xs244nm 

0,(0  =  0) 

[P6A] 

Atkinson 

N02+hv  — N0+0(3P) 

k<398nm 

7.42E-4  s'1 

[P6B] 

Atkinson 

H20+hv-H2+0(3p) 

X<243nm 

0,(0  =  0) 

[P7A] 

Atkinson 

H20+hv-H*40H* 

X<239nm 

1.09E-5  s"1 

[P7B] 

Atkinson 

H20+hv-H2+0(lD) 

X<176nm 

0 

[P7C] 

Atkinson 

C02+hv-CO+0(lD) 

negligible 

C02+hv— C0+0(3P) 

negligible 

HONO+hv  —  NO+OH* 

>w<591nm 

2.61E-3  s'1 

[P50A] 

Atkinson 

HONCM-hv  -*H*+N02 

>.<367nm 

0,(0  =  0) 

[P50B] 

Atkinson 

HONO+hv  — HNO*+0(3P) 

?,<283nm 

0,(0  =  0) 

[P50C] 

Atkinson 

Radical  Chemistry: 
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0H*4<>3=H02*4O2 

6.7*10- 

14 

cm3  molecule- 

1  s" 

1         [Rl] 

DeMore 

H02*+03=OH*+202 

2.0*10" 

IS 

cm3  molecule- 

*  s" 

1         [R2] 

DeMore 

NO+O3=N024O2 

1.8*10" 

14 

cm3  molecule- 

l  s- 

1         [R3] 

DeMore 

NQ24O3=N03-fO2 

1.2*10- 

13 

cm3  molecule"  ■ 

'  s" 

L         [R4] 

DeMore 

HC0*4O2=HQ2*4<X) 

5.6*10" 

12 

cm3  molecule" • 

L  $■• 

L            [R9] 

Anatasi 

HCHO0H*=H*-fCCM-H20 

1.1*10- 

11 

cm3  molecule"- 

'  s" 

1          [R10] 

Seinfeld 

OH*+OH*=H200(1D) 

1.9*10" 

12 

cm3  molecule'- 

l  s" 

L         [Rll] 

DeMore 

0H*+H02*=H2O+O2 

1.1*10- 

10 

cm3  molecule- 

1  s* 

1         [R12] 

Atkinson 

OH*+H202=H2OH02* 

1.7*10" 

[2 

cm3  molecule"- 

•  s" 

1         [R13] 

DeMore 

N02+OH*-HNQ3 

1.1*10-' 

11 

cm3  molecule" - 

■  s"- 

l          [R14] 

Seinfeld 

CO+OH*=C02+H* 

2.9*10" 

13 

cm3  molecule"- 

l  s" 

L         [R15] 

DeMore 

N+OH*=NOH* 

4.9*10" 

11 

cm3  molecule"- 

•  s" 

1          [R16] 

Atkinson 

HCHO0H*=H2OHG0 

1.0*10- 

11 

cm3  molecule" 

1  s" 

1         [R17] 

Atkinson 

H*+02  -  HO2* 

1.2*10"' 

[2 

cm3  molecule" j 

L  s"^ 

L           [R18] 

Anatasi 

H*+03  -  OH*  +  O2 

2.9*10"' 

11 

cm3  molecule" 

1  s" 

I          [R19] 

Anatasi 

H02*+HQ2*=H2Q2+Q2 

1.6*10"] 

2 

cm3  molecule"^ 

■  s"^ 

L           [R20] 

DeMore 

NOH02*=N02+OH* 

8.3*10-] 

L2 

cm3  molecule"^ 

L  s"^ 

L           [R22] 

DeMore 

N02+H02*=H0NO+O2 

1.4*10"] 

14 

cm3  molecule"- 

L  s"^ 

L           [R23] 

Stockwell 

CC>H02*=C02^OH* 

1.9*10** 

52 

cm3  molecule"- 

■  s" 

I           [R24] 

Atkinson 

0(1D)+N2=0(3P)+N2 

1.8*10" 

11 

cm3  molecule"- 

■  s" 

L           [R30] 

DeMore 

0(3p)+02+M=03-fM 

2.9*10'] 

1 

cm3  molecule" - 

•  s~l 

[R31] 

Seinfeld 

NO+0(3p)+M=NQ2->M 

3.0*10"] 

1 

cm3  molecule" ^ 

•  s" 

[R32] 

Atkinson 

NO24O(3p)=02+N0 

9.7*10-] 

2 

cm3  molecule" ^ 

■  s"1 

[R33] 

DeMore 

0(3p)-fO3=024O2 

9.5*10"] 

5 

cm3  molecule" - 

L  s" 

L          [R35] 

DeMore 

0(lD)+H200H*40H* 

2.2*10"3 

0 

cm3  molecule"- 

L  s"- 

I          [R36] 

DeMore 

0(  1D)+03=0+0(3P)+0(3P) 

1.2*10"] 

1 

cm3  molecule"^ 

s"! 

[R37] 

DeMore 

0(lD)+O2=0("3P)+02 

4.0*10"] 

1 

cm3  molecule*- 

^  S" 

L          [R38] 

DeMore 

N+02=NOO(3P) 

8.9*10"] 

.7 

cm3  molecule"^ 

L  S" 

-          [R39] 

Atkinson 

N-fO3=N0+O2 

1.0*10"] 

.6 

cm3  molecule"^ 

^  s"] 

L          [R40] 

Atkinson 

NOOH*=HONO 

i.i*io-] 

1 

cm3  molecule"! 

s"! 

[R51] 

Atkinson 

HON040H*=H2ONCb 

4.9*10"] 

2 

cm3  molecule"  ^ 

s"! 

[R52] 

Atkinson 

CH3OH+OH*=0.85{CH2OH*l+ 

7.8*10"] 

3 

cm3  molecule"! 

s"1 

[R60] 

Atkinson 

0.15{CH3O*+H2O} 

CH20H*+0?=HCHOH02* 

9.8*10-] 

2 

cm3  molecule"^ 

S"! 

[R61] 

Atkinson 

CH30*+02=HCHOH02* 

1.9*10"] 

5 

cm3  molecule"! 

■  s"1 

[R62] 

Atkinson 

Formaldehyde  Reactions   (HCHO): 
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Photolysis: 

HCHO+hv-*H*+HCO* 

X<300nm           1.30E-5S-1 

[P1A] 

Seinfeld 

HCHO+hv—COfH2 

X*300nm            3.86E-5  s-1 

[P1B] 

Seinfeld 

Radical  Chemistrv: 

HCHOOH*=H*+COH20 

1.1*10-11  cm3  molecule-1  s"1 

[RIO] 

Seinfeld 

HCHO0H*=H2OHG0* 

1.0*10-11   cm3  molecule-1  s'1 

[R17] 

Atkinson 

CH2OH*-h02=HCHOH02* 

9.8*  10" 12  cm3  molecule-1  s-1 

[R61] 

Atkinson 

CH30*+02=HCHOH02* 

1.9*10-15  cm3  molecule-1  s-1 

[R62] 

Atkinson 

d[HCHO]/dt=-PlA-PlB-R10-R17+R61+R62 


Methanol  Reactions  (CH3OH): 

Photolysis: 
CH30H-rhv—  negligible 


Radical  Chemistrv: 
CH3OH+OH*=0.85  {CH2OH*}+ 
0.15fCH3(5*+H20i 


7.8*10-13  cm3  molecule-1  s-1 


[P2] 
[R60]      Atkinson 


d[CH3OH]/dt=-R60 


Carbon  Dioxide  Reactions  (CO2): 


Radical  Chemistrv: 
CO  +  OH*-*C02  +  H* 
CO  +  H02*-C02  +  OH' 


2.2*10-13  cm3  molecule-1  s-1 
1.9*10-32  cm3  molecule-1  s-1 


[R15]       DeMore 
[R24]      Atkinson 


d[C02]/dt=P15+P24 


Atomic  Oxygen  Reactions  (0): 
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Photolysis: 


Q3+hv-0(lD)4O2 

Xs411nm 

1.12E-1 

s-1 

[P3A] 

Atkinson 

03+hv-*0(3p)+02 

/,<1180nm 

1.32E-2 

s"l 

[P3B] 

Atkinson 

02+hv-*O(lD)+O(3p) 

X<175nm 

0 

[P4A] 

Atkinson 

02+hv-*O(3p)+O(3p) 

X<242nm 

1.23E-9 

s-1 

[P4B] 

Atkinson 

02+hv-*0(lD)+0(lD) 

X<137nm 

0 

[P4C] 

Atkinson 

N02+hv—NO+O(lD) 

Xs244nm 

0,(O  = 

0) 

[PGA] 

Atkinson 

N02+hv^NO+0(3p) 

X<398nm 

7.42E-4 

s-1 

[P6B] 

Atkinson 

H20+hv  -*H2+0(3P) 

X<243nm 

0,(<t>  = 

0) 

[P7A] 

Atkinson 

H20+hv^H2+0(lD) 

X<176nm 

0 

[P7C] 

Atkinson 

H2Q2+hv-H20+0(lD) 

>.<359nm 

0,(0  = 

0) 

[P5B] 

Atkinson 

HONO+hv  -*HNO*+0(3P) 

X<283nm 

0,(0  = 

0) 

[P50C] 

Atkinson 

Radical  Chemistrv: 

0H*+0H*=H2O+O(lD) 

1.9*10"12  cm3 

molecule"*  s" 

.  a 

[Rll] 

DeMore 

0(1D)+N2=0(3P)+N2 

1.8*10-11  cm3 

molecule"^  s" 

[R30] 

DeMore 

0(3p)+02+M=O34M 

2.9*10-n   cm3 

molecule"!  s"- 

[R31] 

Seinfeld 

NOO(3p)+M=NOz+M 

3.0*10"u  cm3 

molecule"!  s~ 

[R32] 

Atkinson 

N02+0(3P)=02+NO 

9.7*10"12  cm3 

molecule"!  s"^ 

[R33] 

DeMore 

0(3p)+03=202 

9.5*10"15  cm3 

molecule"!  s"^ 

[R35] 

DeMore 

0(lD)+H2O=20H* 

2.2*10-10  cm3 

molecule"!  s"- 

[R36] 

DeMore 

0(lD)4O3=02+20(3p) 

1.2*10-1]   cm3  ] 

molecule"!  s"J 

[R37] 

DeMore 

0(lD)+O2=0(3P)+O2 

4.0*10"1  l   cm3  i 

molecule"!  s"l 

[R38] 

DeMore 

Nht02=NOO(3P) 

8.9*10"17  cm3 

molecule"!  s" 

[R39] 

Atkinson 

aAssumed  the  products  of  Rll  were  negligible. 
Use  Steady  State  Assumption: 

[O(3p)]ss=(P3B+P4A-rP4B+P6B+P7A+P50C+R30+R37+R38+R39)/(KR331 
N02*KR35*03+KR32*NO+KR3 1*02) 

[O(1D)]Ss=(P3A+P4A+P4C+P6A+P7C+P5B+Rll)/(KR30*N2+KR36*H2CM- 
KR37*03+KR38*02) 


NO  Reactions: 
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Photolysis: 
N02+hv-*N0+O(lD) 
N02+hv-*NO+0(3p) 
HONOhv-NO+OH* 

Radical  Chemistry: 
NO  +  O3  =  NO2  +  O2 
N  +  OH*  =  NO  +  H* 
NO  +  H02*  =  N02  +  OH* 
N0  +  0(3P)  +  M=N02  +  M 
N02  +  0(3P)  =  02  +  N0 
N+02=NOO(3P) 
N-f03=NO+02 
NO+OH*=H0N0 


Xs244nm           0,  (<l>  =  0) 

[PGA] 

Atkinson 

\<398nm           7.42E-4  s"1 

[P6B] 

Atkinson 

X<591nm           2.61E-3  s"1 

[P50A] 

Atkinson 

1.8*10"14  cm3  molecule"1  s"1 

[R3] 

DeMore 

4.9*10-11   cm3  molecule"1  s"1 

[R16] 

Atkinson 

8.3*10"12  cm3  molecule'1  s"1 

[R22] 

DeMore 

3.0*10-11   cm3  molecule"1  s"1 

[R32] 

Atkinson 

9.7*10"12  cm3  molecule"1  s"1 

[R33] 

DeMore 

8.9*10"17  cm3  molecule"1  s'1 

[R39] 

Atkinson 

1.0*10"16  cm3  molecule"1  s'1 

[R40] 

Atkinson 

1.1*10-11  cm3  molecule"1  s"1 

[R51] 

Atkinson 

d[NO]/dt=P6A+P6B+P50A-R3+R16-R22-R32+R33+R39+R40-R51 


Hydrogen  Peroxide  Reactions  (H2O2): 


Photolvsis: 

H202  +  hv-»OH*  +  OH* 
H2O2  +  hv-*H20  +  0(1D) 
H202  +  hv— H*+H02* 
H202  +  hv-2H*+02 

Radical  Chemistry: 
0H*  +  H202  =  H?0  +  H0?* 
H02*  +  H02*=H202+O2 


\<577  nm 
X<359nm 
X<324  nm 
X<213  nm 


8.52E-4  s"1 
0,(ct>  =  0) 
0,(0-0) 
0,(0  =  0) 


1.7*10~12  cm3  molecule"1  s'1 
1.6*10"12   cm3  molecule"1  s"1 


[PSA]  Atkinson 

[P5B]  Atkinson 

[P5C]  Atkinson 

[P5D]  Atkinson 


[R13] 
[R20] 


DeMore 
DeMore 


d[H2O2]/dt=-P5A-P5B-P5C-P5D-R13+R20 


NO  2  Reactions: 
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Photolysis: 

N02+hv-*N0+O(lD) 

N02+hv-NOO(3p) 

HONOhv-H*+NQ2 

Radical  Chemistry: 
N0+C^=N024O2 

NQ2+03=N03+02 
N02+OH*=HN03 

NOH02*=N02^0H* 

N02+H02*=H0NO02 

NO+O(3p)+M=N02+M 

NO2+O(3p)=O2+N0 
HONO+OH*=H2ON02 


ks244nm           0,  (<D  = 

0) 

[PGA] 

Atkinson 

X<398nm           7.42E-4 

S"1 

[P6B] 

Atkinson 

>v<367nm           0,  (O  = 

0) 

[P50B] 

Atkinson 

1.8*10"14  cm3  molecule"1  s": 

[R3] 

DeMore 

1.2*10"13   cm3  molecule"1  s" 

[R4] 

DeMore 

1.1*10_1 1   cm3  molecule*1  s"- 

[R14] 

Seinfeld 

8.3*10"12  cm3  molecule'1  s" 

[R22] 

DeMore 

1.4*10"12  cm3  molecule"1  s'^ 

[R23] 

Stockwell 

3.0*1CH  ]   cm3  molecule"1  s" 

[R32] 

Atkinson 

9.7*10"12  cm3  molecule"1  s"- 

[R33] 

DeMore 

4.9*10"12  cm3  molecule"1  s"J 

[R52] 

Atkinson 

d[NO2]/dt=-P6A-P6B+P50B+R3-R4-R14+R22-R23+R32-R33+R52 


Nitrous  Acid  Reactions  (HONO): 


Photolvsis: 

HONOhv-*NO+OH*  " 

X<591nm 

2.61E-3  s"1    [P50A]    Atkinson 

HONO+hv-H*+N02 

>.<367nm 

0,  (<t>  =  0)     [P50B]     Atkinson 

HONO+hv- HNO*+0(3P) 

>.<283nm 

0,  (<j>  =  0)     [P50C]     Atkinson 

Radical  Chemistry: 

N02+H02*=H0NO+O2 

NO+OH*=H0N0 

HONOOH*=H2ON02 


1.4*10"14  cm3  molecule"1  s"1 
l.l*10"n  cm3  molecule"1  s"1 
4.9*10"12  cm3  molecule"1  s"1 


[R23]  Stockwell 
[R51]  Atkinson 
[R52]      Atkinson 


d[HONO]/dt=-P50A-P50B-P50C+R23+R51-R52 


Hydroxyl  Radical  Reactions  (OH*): 
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Photolysis: 
H202+hv-OH*+OH* 
H20+hv-*H*+OH* 
HONOhv-NOOH* 

Radical  Chemistry: 

0H*4O3=H02*+O2 

H02*+O3=0H*+202 

HCHO+OH*=H*4<:C>fH20 

0H*+OH*=H20+0(lD) 

0H*+H02*=H2O02 

0H*+H202=H2OH02* 

N02+OH*=HN03 

COOH*=C02+H* 

N+OH*=NO+H* 

H*+03  =  0H*+02 

NOH02*=N02-MDH* 

COH02*=C02+OH* 

0(lD)+H2O=OH*+0H* 

NO+OH*=HONO 

HONO+OH*=H2ON02 
CH30H+OH*=6.85  { CH2OH*}+ 
0.15{CH36*+H2O} 

b  Assumed  [OH*]  to  be  low,  considered  Rll  to  be  negligible. 

Use  Steady  State  Assumption: 

[OH*]SS=(P5A+P7B+P50A+R2+R19+R22+R24+2'rR36)/(KRl*O3+KR10* 
HCHO+KR12*H02+KR13*H202+KR14*N02+KR15*CO+KR16*N+ 
KR5  l*NO+KR52*HONO+KR60*CH3OH) 


X<557nm           8.52E-4  s"1 

[P5A] 

Atkinson 

X<239nm           1.09E-5  s"1 

[P7B] 

Atkinson 

X<591nm           2.61E-3  s"1 

[P50A] 

Atkinson 

6.7*10"14 

cm3  molecule"- 

'  s"1 

[Rl] 

DeMore 

2.0*10-15 

cm3  molecule"- 

l  s"1 

[R2] 

DeMore 

Lino-11 

cm3  molecule"^ 

■  s"1 

[RIO] 

Seinfeld 

1.9*10"12 

cm3  molecule"- 

s-lb 

[Rll] 

DeMore 

Lino-10 

cm3  molecule"  ^ 

'  s"1 

[R12] 

Atkinson 

1.7*10"12 

cm3  molecule"^ 

s"1 

[R13] 

DeMore 

l.mo-11 

cm3  molecule"^ 

■  s"1 

[R14] 

Seinfeld 

2.2*10"13 

cm3  molecule"! 

•  s"1 

[R15] 

DeMore 

4.9*10-1  ] 

cm3  molecule"J 

■  s"1 

[R16] 

Atkinson 

2.9*10"n 

cm3  molecule"! 

•  s"1 

[R19] 

Anatasi 

8.3*10"12 

cm3  molecule"! 

■  s"1 

[R22] 

DeMore 

1.9*10"32 

cm3  molecule"- 

I-  s"1 

[R24] 

Atkinson 

2.2*10"10 

cm3  molecule"! 

s"1 

[R36] 

DeMore 

1.1*10"12 

cm3  molecule'^ 

s"1 

[R51] 

Atkinson 

6.6*10"12 

cm3  molecule"^ 

s'1 

[R52] 

Atkinson 

7.8*10"13 

cm3  molecule"1 

s'1 

[R60] 

Atkinson 

Hydro peroxyl  Radical  Reactions   (HO2*): 
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Photolvsis: 

H202+hv-*H*+H02* 

X<324nm           0,  (0=0)        [P5C] 

Atkinson 

Radical  Chemistrv: 

0H*-fO3=H02*+O2 

6.7*10"14  cm3  molecule"1  s': 

L         [Rl] 

DeMore 

H02*+03=OH*+202 

2.0*10-15  cm3  molecule"1  s"] 

[R2] 

DeMore 

HCO*+02=HQ2*+CO 

5.6*10*12  cm3  molecule"1  s": 

[R9] 

Anatasi 

OH*+H202=H2CM-HQ2* 

1.7*10"12  cm3  molecule"1  s"] 

[R13] 

Seinfeld 

H*+02=H02* 

1.2*10~12  cm3  molecule*1  s*1 

[R18] 

Anatasi 

H02*+H02*=H202+Q2 

1.6*10~12   cm3  molecule"1  s"- 

l  c       [R20] 

Seinfeld 

NO+H02*=N02+OH* 

8.3*10"12   cm3  molecule*1  s": 

L          [R22] 

Seinfeld 

N02+H02*=H0NO+O2 

1.4*10*14  cm3  molecule*1  s": 

1          [R23] 

Stockwell 

C0+H02*=C024OH* 

1.9*10*32  cm3  molecule"1  s": 

L           [R24] 

Atkinson 

cAssumed  [HO2]  to  be  low,  considered  R20  to  be  negligible. 

Use  Steady  State  Assumption: 

[H02*]ss=(P5C+Rl+R9+R13+R18)/(KR2*03+KR22*NO+KR23*N02^KR24*CO) 


Carbon  Monoxide  Reactions  (CO): 


Photolvsis: 
HCHOhv-CO+H2 

Radical  Chemistrv: 
HC0*+O?=H02*+C0 
HCHO+OH*=  H*+COH20 
COOH*=C02+H* 
CCH-H02*=C02+OH* 


X>300nm 


3.86E-5  s*1      [P1B]     Atkinson 


5.6*10~12  cm3  molecule"1  s"1 

l.l*10"n  cm3  molecule"1  s"1 

2.2*10'13  cm3  molecule*1  s"1 

1.9*10"32  cm3  molecule"1  s"1 


fR9]  Anatasi 

[RIO]  Seinfeld 

[R15]  DeMore 

[R24]  Atkinson 


d[CO]/dt=PlB+R9+R10-R15-R24 


Atomic  Hydrogen  Reactions  (H*): 
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Photolysis: 

HCHOhv-H*+HC0* 

X<300nm           l.30E-5s  -1 

[P1A] 

Atkinson 

H202+hv-H*+H02* 

X<324nm           0,  (<J>  =  0) 

[P5C] 

Atkinson 

H202+hv-2H*+02 

>.<213nm           0,(0  =  0) 

[PSD] 

Atkinson 

H20+hv  -H*+OH* 

X<239nm           1.09E-5  s"1 

[P7B] 

Atkinson 

HONO+hv— H*+N02 

X<367nm           0,  (<J>  =  0) 

[P50B] 

Atkinson 

Radical  Chemistry: 

HCHO0H*=H*+<:OH20 

l.l1 

*10-11   cm^  molecule"1  s"1 

[R10] 

Seinfeld 

CO+OH*=C02+H* 

2.91 

lO"1^   cm^  molecule"1  s"1 

[His] 

DeMore 

N+OH*=NO+H* 

4.91 

"10"11   cm^  molecule"1  s"1 

[rig] 

Atkinson 

H*+02  =  HO2* 

1.2' 

*10"12   cm^  molecule"1  s"1 

[Ri8] 

Anatasi 

H*+03  =  OH*+  O2 

2.9' 

"10"11   cm^  molecule"1  s"1 

[R19] 

Anatasi 

d[H*]/dt=PlA+P5C+P5B+P7B+P50B+R10+R15+R16-R18-R19 


Formyl  Radical  Reactions  (HCO*): 


Photolysis: 
HCHOhv-H*+HC0* 


>,<300nm 


l.30E-5s _1      [P1A]     Atkinson 


Radical  Chemistry: 
HCCr+02  -  HO2*+O0 
HCHOOH*=  H2OHCO* 


5.6*10"12   cmJ  molecu  -*"1  s"1 
1.0*10"n    cm-3  molecule"1  s'1 


[R9]        Anatasi 
[R17]      Atkinson 


d[HCO*]/dt=PlA-R9+R17 


Atmosphe 

ric  Gases:  [Seinfeld,  p.  8] 

Gas 

Ave.  Concentration fDDml 

Ar 

9340 

Ne 

18 

Kr 

1.1 

Xe 

0.09 

N2 

780,840 

Oz 

209,460 

CH4 

1.65 

OQ2 

332 

CD 

0.05-0.2 

H7 

0.58 

N2O 

0.33 

S02 

10-5-10-4 

NH3 

10-4-10-3 

NONO2 

10-6-10-2 

O3 

10-2-10-1 

HNO3 

10-5-10-3 

CH 

0.406*10-5  ppb 

HO2 

0.406*10-2  ppb 

Q3 

30  ppb 

H?0? 

10  ppb 

H?0~ 

22,550 

He 

5.2 
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